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1 Introduction 

CO, 

Lattice-QCD theory predicts the disappearance of the hadronic phase of matter 
once the energy density exceeds a critical value of about 1 GeV per fm 3 [1], giving 
rise to a continuous, deconfined QCD state that is governed by the elementary 
Q . interaction of quarks and gluons. To recreate this phase in the laboratory one 

Q | collides heavy nuclei at relativistic energy with the goal of ascertaining the QCD 

predictions, and to pin down the decay point from the partonic to the hadronic 
phase by obtaining estimates for the transition temperature and energy density. 
The CERN SPS Lead ( 208 Pb) beam facility provides for a top energy of 158 GeV 
per projectile nucleon, corresponding to a total internal CM-system energy of about 
CO ■ 3.5 TeV, to heat and compress the primordial reaction volume. A tenfold higher ^/s 

is reached at the BNL heavy ion collider RHIC. In fact calorimetric data [2] show 
that the average transverse energy density exceeds about 2.5 GeV/fm 3 already at 
CO ■ the SPS in central Pb+Pb collisions. Moreover the study of J/^f production [3] 

demonstrates a suppression of the yield in such collisions, characteristic of the QCD 
"Debye" screening mechanism expected in a deconfined partonic medium [4]. 

This lecture will deal with physics observables that could tell us about that 
medium. As an introduction I will look at the reaction dynamics in AA central 
(— | \ collisions and then turn to my main topic, bulk hadron production and the location 

of the QCD phase boundary. We shall see that the production rates of all strange 
hadronic species play a prominent role, as a signal that helps to understand dense 
QCD matter, and its decay. 
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1.1 Fireball Matter Dynamics 

The initial Woods-Saxon nucleon density distributions, with average energy den- 
sity of about 0.16 GeV per cubic fermi, impinge onto each other leading, at first, 
to nucleon-nucleon collisions occuring concurrently at the microscopic level as en- 
veloped by the overall impact geometry of the target /projectile nuclei. This first 
generation of binary nucleon-nucleon encounters will involve essentially all the inci- 
dent 2 A nucleons if we consider head-on collisions of nuclei with mass number A. In 
Pb+Pb collisions we thus have about 200 primary such nucleon-nucleon collisions of 
first generation. If we could now somehow stop the reaction dynamics, letting the 
reaction products escape to a detector system that identifies them, we would expect 
to register a final multiparticle state with a composition closely resembling A-times 
the (well known) outcome of a nucleon-nucleon collision at similar center-of-mass 
energy ^/s: nothing new! 

However what makes relativistic nuclear collisions interesting is the fact that, in 
reality, the first generation set of A microscopic nucleon-nucleon interaction systems 
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will immediately re-interact while still in statu nascendi of its "pending" output of 
asymptotically distinguishable, "on shell" reaction hadrons which would consist of 
about 10-20 produced mesons and baryon-antibaryon pairs. However, before being 
fully formed they run into a secondary generation cycle of subsequent collisions 
within the nuclear density distributions of the heavy nuclei. In fact there may be 
up to 6 secondary, subsequent collision generations of such pre-exited (not finally 
formed) microscopic collision volumes in a central Pb+Pb collision. Moreover they 
will occur in successive time steps spaced by fractions of a fm/c only - due to the 
relativistic Lorenz contraction of the nuclear density profiles of the Pb nuclei in 
beam direction (7cm = 9 at the SPS). 

Something new! As we have no experimental knowledge of such secondary colli- 
sions of partly incompletely formed ("off shell"), partly decomposed hadrons we have 
to withstand the temptation to capture the overall reaction dynamics in a classical 
billiard ball cascade of subsequent generations of an inelastically multiplying gas of 
known hadrons. Into such a model we would insert the known, vacuum elementary 
cross sections at each microscopic binary encounter and proceed via Monte Carlo 
probability sampling methods. Such approaches are called microscopic hadron trans- 
port models [5]. From the above we would be surprised if they could give a correct 
description of the final outcome of a highly relativistic heavy ion collision because 
their set of microscopic degrees of freedom (isolated binary collisions of "on shell" 
hadrons) fails to capture the unknown nature of secondary, tertiary etc. encounters 
in the dense medium. These encounters, moreover, occur unresolved in time which 
should lead to quantum mechanical coherence, which might render the entire picture 
of isolated, sequential microscopic cascade-sub-processes obsolete. Let us, therefore, 
conclude that the overall large interaction volume of a central Pb+Pb collision will 
(after about 2-3 fm/c of interpenetration and reaction time have elapsed) be com- 
posed of a hitherto unknown state of strongly interacting matter that, however, 
contains all the quantum numbers and relative center of mass energy of the initial 
nuclear projectiles. We thus suspect that it will, quite generally, feature a high spa- 
tial density and, similarly, a high energy density. We therefore sometimes call this 
short-lived object a "fireball". The state of strongly interacting matter inside it is 
the topic of modern QCD theory [6, 7] and of our experiments. 

Now there is something simple about this picture of dense matter in a fireball 
that experiments can check quickly. A new "state" of "matter" supposedly has 
been formed by fusing spheres of target and projectile cold nuclear matter which 
are initially located at opposite ends of longitudinal momentum space in the CM- 
system. By the symmetry of Pb+Pb collisions, a fused fireball should occupy a 
common momentum space volume centered in the CM-system at zero longitudinal 
and transverse momentum. Fig.l shows the distribution in longitudinal phase space 
(measured here by the rapidity variable y = 0.5Zn[(l + (3l)/(1 — (3l)] for the negative 
and neutral K-mesons, and for the anti-strange hyperon A [8]: all are Gaussians well 
centered at ycu = 0. Note that the K~ and A consist only of newly created quarks 
that were not brought into the fireball as initial nucleon valence quarks. For an 
ideal, spherically symmetric fireball we would also require an isotropic momentum 
distribution. Looking simultaneously at the transverse momentum (not shown here) 
and at the y-distributions of these particles we see that this is not strictly the case 
(the fireball is longitudinally stretched and looks like a fire-football of US vintage) - 
but closely enough [9]. 
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Figure 1: Rapidity distributions of K , K and A produced in central Pb+Pb 
collisions at top SPS energy (158 GeV per nucleon). From Ref.[8]. 



Furthermore we can experimentally gather almost all hadrons created in a central 
Pb+Pb collision at top SPS energy: it is about 2500 of them! Superimposing the 
total transverse energy in the football, carried by the various hadronic species we 
can estimate [10] its total average transverse energy density in space. This is the 
2.5 GeV/fm 3 that we mentioned initially [2]. The fireball should thus initially exceed 
the QCD phase boundary between hadrons and quark-gluon matter, which is in the 
vicinity of 1 GeV/fm 3 ! Then it will expand explosively, back to hadrons via the 
phase transformation. 

1.2 Deconfinement and Phase Boundary: Signals 

How do we get experimental signals that elucidate the state of matter in the fire- 
ball maximum density stage? Like with supernova analysis, we may try two ways: 
observe the primordial "light curve" (i.e. electromagnetic and neutrino radiation) 
or look at the bulk explosion material (e.g. expansion modes, element composition 
etc.). For the nuclear fireball the neutrino signal is out but we are left with "black 
body radiation" : directly emitted photons and lepton pairs which leave the fireball 
while it is hot and dense. The corresponding experiments [11] are very demanding 
but successful. However I will not cover them here for lack of space. 

Focusing on experiments that analyze the fireball material after expansion and 
cooling I briefly mention the suppression of the J/^f yield confirmed experimentally 
[3] to occur in central Pb+Pb collisions at top SPS energy, 158 GeV per nucleon. 
Still discussed vigorously in the community, we may take such observations to be, at 
least, consistent with creation of a fireball energy density well in excess of 1 GeV/fm 3 
at which QCD predicts deconfinement thus dissolving the J/^f. This observation 
agrees with our above estimate of the transverse energy density as derived from the 
total final phase space density of the produced hadrons. 

If one tentatively takes for granted such indications of a deconfinement state 
at top SPS energy one expects, likewise, to receive signals of the parton to hadron 
phase transformation bound to occur once the primordial high density state expands 
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and cools toward the critical temperature and energy density. This brings me to the 
principal topic of the following chapters, in which I shall develop this observable in 
broader detail: could it be imagined that the detailed composition of the fireball 
decay products, in terms of abundances of the various hadronic species, captures 
the parton to hadron QCD phase transformation period of the dynamical fireball 
evolution? Or, more generally speaking: we expect to create, at least, a novel state 
of strongly interacting matter in the short lived fireball which must then decay We 
infer from the above considerations of microscopic fireball dynamics that this high 
density state should feature novel degrees of freedom at the microscopic level. These 
may be partons (at top SPS energy and beyond), or shadows of the familiar hadron 
spectrum as modified by interaction with their dense, surrounding fireball medium 
(as might occur at lower collisions energies), and by the microscopic reaction cycles 
during the interpenetration stage of the nuclear density distributions. In any case it 
is the decay of this unknown coherent quantum mechanical state to a quasi classical 
state with familiar degrees of freedom, in our case a hadron and resonance gas, that 
we have to consider. The composition of that hadron gas is detectable: our signal. 

Let me make a few remarks referring to QCD folklore. Almost trivially, all con- 
ceivable states of strongly interacting matter are falling under QCD governance. 
Thus in a decay of a state higher in energy density and temperature, upon fireball 
expansion, to the lower density hadron gas it is just the set of proper QCD degrees 
of freedom that is changing - we are thus dealing with a QCD phase transformation. 
At the hadron gas level, QCD resides in the spectrum of mass, spin, isospin, flavour, 
width etc. of the vast array of hadronic species. Thus in a QCD parton to hadron 
phase transformation, an initial fireball ensemble of flavoured plus coloured quarks 
and antiquarks, and of colour-anticolour carrying gluons is, sloppily speaking, "look- 
ing down" at the QCD realization below: the hadron spectrum. The ensuing, colour 
neutralizing "condensation" process of partons leads to hadronic objects featuring 
a spectrum of specific flavour-colour-spin-momentum internal compositions that ab- 
sorb the initial partonic degrees of freedom. Physics experience suggests that such 
a process, occuring in an extended volume, should be governed by statistics. In the 
partonic view "from above" a light (ud positive) pion will be certainly the easiest 
way to hadronize in a predominant u, u, d, d population (likewise for the other 
pions). In strong contrast, a heavy {sss) Omega hyperon will be a highly unlikely 
hadronization outcome, in view of competing K~ and K° mesons that could take 
care of the s-quarks combining them with the more abundant u and d quarks. Net 
result of all of this: the hadronization transition populates the hadronic spectrum 
in order of the relative statistical weight of the hadrons. We call such a process 
"phase space (statistical weight) dominated". It will thus create a hadron gas en- 
semble of maximum entropy, a decoherent, classical system [12] which exhibits a 
characteristic ordering pattern concerning the relative abundance of each hadronic 
species. It is, thus, not surprising that the abundance spectrum (expressed in terms 
of the various average hadronic multiplicities per collision event) obeys a statistical 
Gibbs ensemble [13], either in the so-called "canonical" or in the "grand-canonical" 
realization. 

The statistical analysis of the composition of multihadronic final states was pio- 
neered by Hagedorn in the 1960's [14]. It was revived in the last decade, applying it 
to the hadron populations observed in relativistic nucleus-nucleus collisions as well 
as in pp and e + e~ annihilation [13]. In this approach one captures the temperature 
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and energy density prevailing at birth of the multihadronic final state, i.e. the point 
in the fireball dynamics where it decouples, by decoherence, from the novel state of 
high energy density/temperature created in the early phase of the dynamics. This 
is clearly an interesting signal! 



2 Hadron Multiplicity and Strangeness Enhancement 

The first SPS experiments with 32 S-beams at 200 GeV/A showed an enhancement 
of various strange particle multiplicities, chiefly K + , A and A, relative to pion mul- 
tiplicities, in going from peripheral to central S + (S, Ag, Au) collisions [15]. This 
observation appeared to be in-line with the pioneering analysis of Rafelski and Miiller 
[16] who first linked strangeness enhancement to the advent of transition from the 
hadronic to a partonic phase. This offered lower effective ss threshold, shorter dy- 
namical relaxation time toward flavour equilibrium, and an additional, nontrivial 
effect of relatively high net baryon number or baryochemical potential: the light 
quark Fermi energy levels move up, perhaps even to the s-quark mass at high 
and the Boltzmann penalty factor for the higher mass ss pair creation might be 
removed. This latter aspect was mostly ignored in the late 1980s but receives fresh 
significance as we become increasingly aware of the crucial role of \xb- 

In this section, and in sections 3 and 4, I will present a sketch of our recent 
progress, both in gathering far superior data and in the understanding of the statisti- 
cal model that was rudimentarily anticipated in such early strangeness enhancement 
speculations. 
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Figure 2: Phi meson to pion multiplicity ratios for Pb+Pb at 158 GeV/A, as a 
function of collision centrality given by the number of participating nucleons (left), 
and by R-b/2 where R is the nuclear radius. Data for min. bias pp and for central 
light nucleus collisions are also given, Ref.17 and 20. 
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Fig. 2 shows a modern version of the typical strangeness enhancement phenom- 
ena. </>-meson to positive pion multiplicity ratios obtained by NA49 [17, 20] in 
Pb+Pb SPS collisions at 158 GeV/A (corresponding to ^=17.3 GeV) are plotted 
for a sequence of collision centrality conditions from peripheral to central. At the 
peripheral end the minimum bias p+p point [18] matches with the trend. The cen- 
trality scale employed here is, at first, the number of participating nucleons (left 
hand side of Fig.2). The raw data centrality bins are ordered in NA49 data by de- 
creasing projectile spectator energy as recorded in a zero degree calorimeter. This 
information is converted to mean participant nucleon number by a VENUS calcu- 
lation from which one also derives the mean impact parameter b, corresponding to 
each centrality bin. Neither of these scales turn out to be satisfactory [17, 19, 20] 
in merging data from central light nuclei collisions such as C+C with the various 
centralities of the Pb+Pb collisions. For example, a central C+C collision has b ~ 2 
fm and N part « 18 but on a b scale the 4>/ir + value is about 50% lower than the 6 = 2 
result for the much heavier Pb+Pb system. Inversely on the N part scale: N part =l8 
corresponds to very peripheral Pb+Pb and the central C+C result is about 50% 
higher than the Pb+Pb curve. A central collision of a relatively light nuclear pair 
thus behaves quite differently from a very peripheral heavy nuclear collision where 
only the dilute Woods-Saxon density tails interact! The scale of the right hand side 
of Fig.2 is an intuitive guess [19] to represent the relative compactness, or volume-to- 
surface ratio of the primordial interaction volume, by the variable R-6/2, where R is 
the radius of the colliding nuclear species. It might be connected with the energy or 
collision density reached in the primordial collision volume (see chapter 5). We see 
that the central light nuclear collision data now merge with the Pb+Pb centrality 
scale. 

The "strangeness enhancement factor" is given oftentimes as the production 
ratio of AA central/(pp min. bias times 0.5 N part ). In the case of Fig.2 it would 
be roughly 2.5. We also have a systematic study of strangeness enhancement in 
Au+Au collisions at the BNL AGS energy of 11 GeV per nucleon where this factor 
is about three [21]. Multistrange hyperons [22, 23] show factors between 4 and 15 
at top SPS energy. 

Bulk strangeness enhancement in central collisions is a nuclear feature, absent in 
pp collisions. Of course we lack a detailed picture about "centrality" in pp collisions 
but we could still employ e.g. the total charged particle multiplicity to select more 
or less "violent" collisions. Fig. 3 shows the K + /ir ratio of pp at 158 GeV versus 
charged particle multiplicity to be essentially flat [24]. Similar findings are made up 
to Tevatron yfs = 1.8 TeV pp collisions: the K + /tt ratio is 50% higher here but also 
almost independent of N c h [25]. 

A picture emerges in which strangeness enhancement, or more generally speaking 
the yield order in the overall bulk hadron population is connected with "sequential- 
ity" of interactions at the microscopic level, i.e. with the fact that a number of 
successive collisions occurs in narrow intervals of time if one may employ a naive 
Glauber picture. This number reflects the size and density of the primordial in- 
teraction zone. It thus seems that a certain (extended) size and a certain inter- 
action rate (and thus energy density), typical of central nuclear collisions, conspire 
in strangeness enhancement. Unfortunately we do not really know how to describe 
a second, third etc. collision of a hadron, within fractions of a fm/c space-time 
distance. We would then have probably resolved the key issue: does it dissolve into 



6 



A 


0.12 






V 




A 


0.1 


u 




V 




0.08 




0.06 




0.04 




0.02 








NA49 preliminary 
i i i I i i i I i i i I i 



2 4 6 8 10 12 14 

N 

ch,measured 



Figure 3: The multiplicity ratio of positive kaons to pions in pp collisions at 158 
GeV, as a function of charged particle multiplicity, Ref.24. 



a parton cascade from which the final hadrons are reconstituted? Proton-nucleus 
collisions should hold a key to this question but nobody has succeeded in isolating 
the second, third, n'th successive collision of the projectile, as of yet [26, 27, 28]. 

At the moment we thus forgo pA as an intermediate step although it certainly 
also features changes in the hadronic production ratios [26, 27, 28] and base the 
analysis on comparing pp to AA. Fig. 4 shows the hadronic multiplicities, from pion 
to cascade hyperon, obtained by NA49 for min. bias pp at v / s=17.3 GeV [29]. The 
data are confronted with the Hagedorn statistical model in its canonical Gibbs en- 
semble form as employed by F. Becattini [30], leading to very good agreement (as 
it is well known also for other elementary collisions and energies [31]). The three 
parameters are T=186±7 MeV, a reaction volume of about 7 fm 3 , and a total of 
about 0.5 ss pairs. The apparent validity of a statistical weight-dominated pic- 
ture of phase-space filling has been considered a puzzle already since Hagedorn's 
time. It is clear, however, that the apparent canonical "hadrochemical equilibrium" 
pattern can not result from " rescattering" of produced hadrons: there is none. In 
Hagedorn's view [32] a creation "from above" must hold the key to the apparent 
maximum entropy state, i.e. the QCD process of hadronization [33]. In it a medium 
of unknown composition (the "fireball" or "string" or "cluster") undergoes a quan- 
tum mechanical decay which is dominated by the phase space weight distribution 
that corresponds to the spectrum of hadrons and resonances [12, 34]. The canonical 
abundance pattern and the T-value are a fingerprint of QCD hadronization - do AA 
collision data at high yfs also confirm this picture (they must, of course, also result 
from a hadronization process)? 
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Figure 4: NA49 data for hadron multiplicities in pp collisions at 158 GeV confronted 
with the canonical model of Becattini, Ref.30. 



3 AA collisions in the Grand Canonical Model 

Fig. 5 shows the grand canonical fit by Becattini to the NA49 data from central 
Pb+Pb at 158 GeV/ A [30]. The temperature is 160 ± 5 MeV and /j, B =240 MeV; 
besides, this model employs the much discussed strangeness undersaturation factor 
7 S =0.8. Leaving the second order concern about 7 S to the theoretical community I 
note here that Braun-Munzinger et al. [35] fit a set of data at the same top SPS 
energy without introducing a 7 S ; they report T=170 ±5 MeV, at /i£=270 MeV, 
close enough. There are also studies of the new RHIC STAR data [36] at y/s=130 
GeV by this model [37] and by Kaneta and Xu [38], averaging at 175 ± 5 MeV 
and // B =48 MeV. And the new data of NA49 [39, 40] at 80 and 40 GeV/ A have 
resulted in fit values of T=155 MeV, /i B =270 MeV and T=150 MeV, /i B =395 MeV, 
respectively [41]. I will return shortly to a further discussion of the grand canonical 
approach but wish to, first of all, show an overall impression from these analyses 
which are confronted in Fig. 6 with the new lattice QCD calculations at finite /is by 
Fodor and Katz [42]. 

The latter predict the T, /i dependence of the QCD phase transformation which 
in this model consists of a crossover for all /J, < 650 MeV, i.e. in the SPS-RHIC 
domain. Note that physics observables can change rapidly in a crossover, too: the 
familiar steep rise of e.g. lattice e/T 4 at T c does not, by itself, reveal the order of the 
phase transformation [6]. Anyhow: the hadronization points from grand canonical 
ensemble analysis merge with the phase transformation site of lattice QCD at top 
SPS and RHIC energy. Quite a remarkable result, but also a plausible one [33, 43] 
if we recall that Ellis and Geiger did already point out in 1996 that hadronic phase 
space weight dominance appears to result from the colour-flavour-spin-momentum 
"coalescence" of partons that occurs at hadronization [12, 44]. Unfortunately a 
rigorous QCD treatment of the parton to hadron transition is still missing. 
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Figure 5: Hadron multiplicities for central Pb+Pb collisions at 158 GeV/A from 
NA49 confronted with the grand canonical statistical model, Ref.30. 



At this point the following objection is always raised: if the same basical model 
describes hadronic yield ratios in pp, e + e~ and in central AA collisions, Figs. 4 
and 5, what is special about AA, as you will not tell us now that a QGP is also 
formed in pp?! Answer: on the one hand both collision systems reveal the QCD 
hadronization process which features, furthermore, the Hagedorn limiting hadronic 
temperature Tjj- At top SPS and at RHIC energy T (hadronic ensemble) « Th ~ 
T c (QCD), this is the common feature; it should not be a coincidence. On the 
other hand hadronization appears to occur under distincly different conditions in 
A A collisions, as captured in the transition from a canonical to a grand canonical 
description. Inspection of Fig. 4 and 5 shows that the hadronic population ratios 
are quite different: the falloff from pions to strangeness- two cascade hyperons in 
the former case is about four orders of magnitude whereas it reduces to three in 
the grand canonical situation: strangeness enhancement! In the canonical case the 
small reaction fireball volume is strongly constrained by local conservation of baryon 
number, strangeness neutrality and isospin whereas these constraints fade away in 
the grand canonical ensemble which represents a situation in which, remarkably, 
these conservation laws act only on the average, over a rather large volume, as 
captured by a collective chemical potential [is- This leaves one global quantity fis 
essentially in charge of all the conservation tasks. Note that the statistical model 
does predict nothing, it merely captures this most remarkable feature of the hadron 
gas emerging after hadronic freeze-out. Its observed success implies some kind of 
long range collective behaviour in the hadronizing source, the origin of which is yet 
unknown, but must be specific to central AA collisions. We shall address this topic 
in the following Sections but note, for now, that Strangeness enhancement is 
the fading away of small volume canonical constraints, in the terminology 
of the statistical model [45]. 

This aspect has been recently studied in all detail by Cleymans, Redlich, Tounsi 
and collaborators [45, 46]. Fig. 7 illustrates their results concerning the transition 
from canonical to grand canonical behaviour with increasing number of participants, 
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Figure 6: The lattice QCD phase boundary in the plane of T vs. fiB, Ref.42. 
The hadronization points captured by grand canonical analysis for SPS and RHIC 
energies are also shown. 



i.e. overall "source" size. It is intuitively clear that it should occur, first, in singly 
strange hadrons, the increase occuring with offset (but having a larger specific effect 
on the yields per participant) in S=2,3 hyperons. 

A further, appropriate critical question: how can we understand the other aspect 
of Fig. 6, i.e. the steep falloff from the QCD transition domain occuring at the lower 
SPS energies? We even have a further GC analysis, at top AGS energy, by Stachel 
[47], for central Si+Au collisions at 14.6 GeV/A, shown in Fig.8. The result is T=125 
MeV, /ie=540 MeV, far below the T scale of Fig. 6. The picture of a direct parton to 
hadron transition is intuitively inapplicable at these lower energies. Still the overall 
dynamical trajectory that ends in hadronic chemical (abundance) freezeout should 
arrive there " from above" as hydrodynamical models [48] show. The hadronic system 
is initially very dense along its trajectory, it is not a hadron gas, and it is thus a 
quantum mechanical coherent state composed of excited and in-medium modified 
hadrons that decays to the finally observed classical hadron gas ensemble. This decay 
might also ignore the classical concept of a relaxation time. Recall the nucleus, also 
still a dense system: we do not invoke relaxation time in a transition within such 
a quantal medium, such as /3-decay. And yet "Fermis Golden Rule" asserts that 
the transition strength depends "only" on the squared matrix element times final 
state phase space volume weight plus global conservation laws. And we know that 
the phase space factor oftentimes far overrides the matrix element, in the net decay 
strength. High density hadronic matter behaviour is essentially unknown: an old and 
new research paradigma but at present we state that it is not a classical hadron- 
resonance gas and that it may decay quantum mechanically, the decay products 
well described by a grand canonical ensemble of hadrons and resonances. At top 
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Figure 7: The canonical to grand-canonical transition as reflected in the canoni- 
cal suppression factor which is the inverse of strangeness enhancement, shown for 
strange hadron species with s=l,2,3 at top SPS energy, Ref.46. 



SPS and at RHIC energy, in turn, the increasingly "explosive" nature of partonic 
and hadronic expansion may almost instantaneously dilute the hadronizing source 
toward chemical freezeout, as indicated by T(GC)« T (Hagedorn) rs T c (QCD). We 
may, thus, directly look here at the QCD parton to hadron phase boundary, located 
at T w 170 MeV. 

4 Working with the Grand Canonical Ensemble 
4.1 The Hagedorn Fireball 

The previous sections have described our recent attempts to understand the order 
of multiplicity per collision in which the various hadronic species are populated. We 
have referred to the statistical Hagedorn model [32] in its modern versions. It has 
turned out that elementary collisions, such as pp, pp, e + e~ — ► hadrons, are well 
described in the canonical version of the statistical model, whereas central nucleus- 
nucleus collision hadronic final states appear to obey grand canonical statistical 
order. Let me emphasize again that in this model we consider the finally observed 
multihadron state to result from the decay of a quantum mechanically coherent 
fireball stage that resides "above" the finally observed, "frozen-out" classical hadron 
gas, in terms of energy density and temperature. Its decay occurs at a certain, late 
stage in the overall dynamical evolution, imbedded into an overall cycle of initial 
interpenetration, fireball stage of maximum energy density, and expansion dilution 
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Figure 8: Hadron yield ratios at top AGS energy, in central Si+Au collisions at 14.8 
GeV/A as fitted with the grand canonical statistical model, Ref.47. 

and cooling. The composition and structure of matter in the intermediate fireball 
stage is the object of foremost interest. Its direct radiation output signals form one 
set of appropriate observable signals [11]. The other signals are derived from the 
fireball freeze-out decay into hadrons. 

The freeze-out product of the fireball is sufficiently dilute to be quantum mechan- 
ically incoherent, thus being describable as a classical hadron gas. The multiplicities, 
and multiplicity ratios of the various hadronic species in this gas form an image of 
the instant of decay: we thus study the decay properties of the unknown state of mat- 
ter in the high density fireball: the conditions prevailing " at birth" of the frozen-out 
hadron gas, common to all hadrons emerging from the fireball. These common con- 
ditions are, chiefly, the temperature, energy density and net baryon number density. 
They are captured by the statistic ensemble analysis. Its success in describing the 
composition of the frozen out hadron gas (residing in the multiplicity ratios of dozens 
of hadronic species, from pions to multiply strange heavy hyperons) shows that the 
entire, at first sight exceedingly complicated final state, of up to several thousands 
of hadrons, can be understood from a few common macroscopic parameters. To 
summarize: it is the underlying assumption of Hagedorn analysis that the "decay 
from above" is dominated primarily by the energy density (and conserved quantum 
numbers) at the instant of decay, and by the statistical weight distribution offered 
by the hadronic spectrum into which the decay occurs. Immediately after decay 
the hadronic gas thus exhibits a statistical equilibrium population, as captured by a 
Gibbs ensemble of mixed hadronic species. At the risk of overdue repetition, I stress 
again that this apparent "thermal" equilibrium is a result of the decay process, the 
nature of which lies well beyond the statistical model which "merely" captures the 
apparent statistical order, prevailing right after decay. The hadrons and their weight 
distribution are "born into" this state [32, 33]! The observed equilibrium is, thus, 
not achieved by inelastic transmutation of the various hadronic species densities, in 
final hadron gas rescattering cascades, i.e. not by hadron rescattering approaching 
a dynamical equilibrium. The system explodes and cools rapidly after the initial 
hadron gas formation phase [33], the canonical or grand canonical order staying 
frozen-in throughout subsequent expansion (while the momentum spectra etc. still 
get modified by elastic interactions, resonance decays etc.). Hadro-chemical compo- 
sition freeze-out thus occurs prior to the final decoupling from all strong interaction 
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(spectral freeze-out). In a typical central Pb+Pb collision at top SPS energy, 158 
GeV per projectile nucleon, hadro-chemical freeze-out occurs at a temperature of 
165 ±5 MeV (thus capturing the hypothetical QCD phase boundary) but final, spec- 
tral decoupling occurs at about 110 MeV [49], as we learn from hadronic spectra 
and correlations. 

I have inserted this somewhat lengthy section to avoid the misunderstandings and 
controversies that have accompanied the Hagedorn model ever since its inception. 
In the 1970's particle physicists were seeking an understanding of proton collisions 
in analogy to the elementary Feynman graphs of QED, thus trying to insolate sim- 
ilarly "elementary" hadronic processes, connecting the initial and final constituent 
(massive) quarks. However the most probable outcome of a pp collision, namely 
to go to a ten-hadron final state (of statistically varying microscopic composition) 
clearly defied such a picture, thus being called "background" reactions, outside of 
the primary research focus. On the other hand Fermi, Landau and Hagedorn just 
put the emphasis on this side of hadron collisions, guided by the intuition that an 
overall process of ever increasing density of potentially coexisting microscopic sub- 
processes should approach the statistical situation: all subprocesses would feed into 
a symbolic, intermediate "compound state" (reminiscent of Bohrs excited nuclear 
compound state)that was called fireball. Its decay would then feed into a frozen-out, 
statistically ordered hadron gas. The appropriate formulation turned out to be a 
canonical Gibbs ensemble [32]. We have seen the modern version of this model in 
Fig.4. 

Furthermore, the Hagedorn concept of freeze-out from the symbolic fireball (of 
microscopically unknown degrees of freedom) to a classical hadron gas of statistical 
weight dominated composition - that then travels on to final observation unchanged 
by the subsequent dynamical evolution - took a long time to be comprehended. On 
the other hand we are cordially familiar with the same picture from explosive Big 
Bang nucleosynthesis. We know that the cosmic average proton to helium compo- 
sition ratio froze out to its observed ratio once the inelastic transmutations in the 
cosmic fireball (among the various light nuclear species) stopped at about 1 MeV 
temperature - thus presenting to us highly relevant data concerning this dynami- 
cal stage. With further expansion this chemical composition travels on unchanged 
(frozen- in), while the spectral temperature of the cosmic inventory has dropped 
down to a few 10~ 4 eV. Likewise, hadro-chemical fireball freeze-out creates hadron 
composition ratios that travel onward throughout further hadronic expansion. 



4.2 The Grand Canonical Ensemble 

Let us now turn to the formalities of grand canonical ensemble analysis. One starts 
from the formulation of a partition function which specifies the relative weight Zi 
(it counts the sum of possible states) for each particle (or resonance) species i in a 
multihadronic mixed gas at temperature T: 

, ? _ 9iV [°° ^dk 

Hl 67r 2 Tio Ei(fc)[ea;p{(Ei(fc)-Mi)/r}±l] U 

where gi is the statistical Lande factor (the "statistical degeneracy") of species 
i, V is the total common volume shared by all species, E? (k) = k 2 + mf the to- 
tal energy square of species i at momentum k, and Hi = fisBi + fx s Si + [ijli the 
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"chemical potential" of species i. The latter is the typical, unfamiliar ingredient of 
the grand canonical ensemble. Please consult textbooks describing its occurence in 
Zj as a result of adding a so called Lagrange multiplier to the Lagrange density of 
the system in order to enforce global conservation of certain net quantum numbers 
that are specific to the total system contained in the total volume V. For our case of 
"hadro chemistry", these are the net quantum numbers initially carried into the fire- 
ball volume by the incident nuclei. Their total baryon number B, total strangeness 
S, and total isospin (Z-component) I are initially given by the total participating 
nucleon (=baryon) number, S being zero (nuclei have zero strangeness), and by the 
net isospin of all participant nucleons: I=(Z-N)/2. These net quantum numbers are 
conserved by strong interaction, and accompany the collision volume throughout the 
fireball evolution, during which they will be re-distributed over the entire ensemble 
of the hadronic species that is being created. Thus the three components of ^ in 
Zj represent the net impact of the overall quantum number conservation on each 
separate species i. We notice, most significantly, that the conservation laws are thus 
enforced on average only, not locally but over the entire fireball volume. This to- 
tal volume thus acts somehow coherently, balancing globally the quantum number 
exchanges that have occured microscopically during fireball dynamics preceding the 
decoupled hadron gas that is described by this set of Zj's. From our considerations 
above, and following Hagedorn's advice [32], we state again that this global coher- 
ence is not achieved by inelastic relaxation toward equilibrium within the frozen-out 
hadron gas ensemble but must be a characteristic feature either of the quantum me- 
chanically coherent high density fireball preceding freeze-out, or by the mechanism 
of its decay process to the decoupled hadron gas (or by both acting together). After 
decay the resulting hadron gas is imprinted with these global coherence aspects of 
its birth process. The grand canonical ensemble description merely "takes note" 
of the freeze-out product of such a process, analyzing a temporary snapshot of the 
apparent order in the frozen-out state, ignorant of its dynamical origin. 

A short comment concerning the terminology: chemical "potential". You see 
from equ. (1) that the exponential "penalty factor" in the denominator, that ex- 
presses the cost of realizing a certain hadron with total energy Ej(k) within a bath 
of temperature T (characterizing the average kinetic energy, or energy density), gets 
modified by /itj. This modification thus takes into account the "affinity" - to employ 
the classical terminology - that the medium offers to species i. In modern terms 
we realize that /ij in the exponent acts like an average potential for species i in the 
medium, modifying its vacuum energy Ei(k), to now read Ei(k) — /ij. The medium 
thus modifies the cost of realizing species i at momentum k. An "affine" medium 
(positive fii) reduces the penalty factor and thus increases the relative weight Zj. 

From the partition function in (1) we derive the distribution of number densities 
of the species i in the medium, 



The contribution of species i to the overall energy density ej is obtained by 
inserting the factor E^k) in the numerator of the integral. From the number density 
of species i we get its total multiplicity per event by multiplying ni with the total 
reaction volume V. 

Now we want to analyze experimental data, with the set of TVj hadronic multi- 





14 



plicities supposed to be known. From this we have to determine the basic parameters 
V, T and hb, fJ-s and Hi- It turns out that we can, before proceeding to a fit to equ. 
(2), reduce the three chemical potentials to one remaining quantity. We do this by 
exploiting global baryon, strangeness and isospin conservation in a set of coupled 
equations which, in the end, express fi s and fj,j in terms of \ig which thus remains as 
the only independent potential, expressing all conservation tasks. Three parameters 
remain: V, T and /j>b- We have shown in Fig. 5 the outcome of such a fit procedure 
as applied to a set of NA49 hadron multiplicity data for central Pb+Pb collisions 
at top SPS energy. Repeating at various energies [41] we get the systematics of 
hadronic chemical freeze-out conditions as represented in the T, \ib plane (Fig.6). 
That is, basically, how it works. 

An easy qualitative understanding of the procedure that determines the two 
parameters T and //g, in a first cycle (which then leads to fixing V in a second step) 
is obtained from considering the special case of particle to antiparticle multiplicity 
ratios. For easy approximation we go from Bose or Fermi statistics (resulting in the 
±1 term in the denominator of equ. (2)) to Boltzmann classical statistics dropping 
that term. We consider the ratio of a certain particle-antiparticle multiplicity, 

< m > I < Hi >= riiV/rliV = rii/rii (3) 

in which the common total volume parameter drops out. Consideration of the ap- 
propriate ratio of two integrals (2) shows that almost everything drops out because 
statistics and phase space are identical for particle-antiparticle pairs. What remains 
is 

< rii > I < Hi >= exp{(m - TIi)/T} = exp(2m/T). (4) 

If we do this for several such ratios in combination, say K + /K~, p/p, A/A, f2/fi 
etc., we get a set of equations from which to determine the maximum likelyhood 
{/Ub, T} combination. In a second step we fix the volume parameter V by fitting 
equ. (2) to the pion multiplicity < > with known //g and T. 

In the real procedure one has to go through one additional complication. The 
finally observed hadrons are not really the ones that existed at freeze-out, T=165 
MeV or so, in the created hadronic gas. It is composed, at first, of excited hadronic 
states such as A*, A, K* etc. and of resonances like g and u. It is their population 
that obeys the set of equations (2)! After being established at freeze-out this pop- 
ulation decays in well known ways to the finally observed ir, K, A, p, p etc.. These 
latter particles have never been in a state that obeys the grand canonical multiplicity 
order! They serve as an observational input to a grand canonical ensemble fit via an 
attached procedure that was invented by Wroblewski [50], which relates the final, 
observed set of hadron multiplicities to a set of excited hadron and resonance multi- 
plicities that forms the initial hadron composition at freeze-out, via the known decay 
branching ratios. This important backward-transformation was oftentimes ignored 
in early applications of the statistical model, causing considerable confusion. 

5 The Onset of Grand Canonical Order in AA Collisions 

We infer from the previous section that grand canonical ensemble analysis captures 
a snapshot at birth of the classical hadron gas, emerging from central collisions 
of heavy nuclei. The apparent most interesting feature was the quantum number 
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conservation on average only, over a large volume of microscopic ingredients. We 
have claimed that this coherence is not a property of the frozen-out hadron gas 
but should reside in its birth process "from above", referring hypothetically, to the 
dynamics of fireball evolution preceding decoupling and freeze-out. We have, thus, 
invoked a prior dynamics: coherence in a large volume, high energy density fireball 
(of unknown matter composition) which, upon decay, creates grand canonical order. 
The grand canonical analysis merely takes note of this order; it does not know about 
our additional, hypothetical speculation concerning its dynamical origin. 

On the other hand, as we are interested, in particular, in just this dynamical ori- 
gin (as it might tell us about the properties of the unknown fireball matter above: its 
decay temperature and chemical potential), we would now want to inspect a wider 
variety of fireball hadronization processes. We thus turn to elementary hadron col- 
lisions, and to e + e _ -annihilation, in which the interaction volume is very small in 
comparison to AA fireball volumes. Nevertheless, hadronization must occur under 
similar QCD dominance here. This leads us to expect that the statistical temper- 
ature, which is the decay image of the critical energy density at which any QCD 
hadronization occurs, must be common to both elementary and AA collisions at 
yfs > 20 GeV (top SPS to RHIC energy). Because we assume, from all existing 
observables, that we are witnessing a parton to hadron QCD phase transformation 
at these high energies. Indeed the temperature parameters all fall within 170 ± 10 
MeV, from e+e~ to Au+ Au [13, 31, 36, 38, 41]. 

The apparent difference between small and large freeze-out volumes appears, 
thus, to arise not from the decay temperature nor from the final hadronic level 
spectrum but from a specific difference in the size of the coherently decaying system. 
Indeed, in the canonical decay mode the "volume" amounts to a few fm 3 only. The 
relevant conserved quantities are forced in the canonical case to be conserved locally 
- if an Omega hyperon is to be produced it must be locally accompanied by an 
anti-Omega or by an anti-Lambda plus two K + and a tt~ (charge, strangeness and 
baryon number conservation), etc.. Of course this requirement is relatively hard to 
meet locally, and it inflicts a high statistical " punishment factor" . We call this effect 
canonical suppression. In the grand canonical case nothing couples to the Omega 
locally, we only require three strange quarks and some energy density. The resulting 
enhancement of the production rate can be understood as an absence of canonical 
suppression. This statement forms the essence of the analysis [46] of the transition 
from C to GC behaviour that we showed in Fig. 7. 

Before looking at the formal expression for this transition, as derived by Tounsi 
and Redlich [46] let me try a crudely oversimplified qualitative argument for illustra- 
tion, involving the punishment factor only. Consider Omega formation by decay of 
a partonic medium. Canonical decay then requires 3 s quarks and 3 s quarks jointly 
participating (plus other factors that I ignore here). The penalty factor arising from 
local strangeness and baryon number conservation will, in this naive picture, be 
exp(—6m s /T) whereas in a GC situation it diminishes to exp(— 3m s /T), the conser- 
vation being taken care of elsewhere in the large volume. The enhancement factor 
then is exp(3m s /T). Taking m s =140 MeV and T=165 MeV we get a factor of about 
13: rather close to the ratio of about 15 for the Omega production per participant 
nucleon pair in central Pb+Pb collisions at top SPS energy, relative to the p+p pro- 
duction rate at similar energy [23]. The Omega is, of course, a very rare species at all 
incident energies. The "bulk" strangeness abundance in the system is more directly 
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reflected in the singly strange K + yield. Let us try: the canonical penalty factor 
would naively be exp(—2m s /T) for a K + , but only exp(—m s /T) in a GC scenario. 
The GC enhancement factor of bulk strangeness would thus be exp(m s /T) = 2.3: 
close enough to the experimentally observed [17] kaon enhancement, by a factor of 
about two. 

Let us take the above as an amusing illustration concerning the Boltzmann 
phase space factor for different hadronic species occuring in the denominator of 
equ. (1). We have, of course, ignored the effect of the baryochemical potential hb', 
and many other appropriate considerations. An exact treatment can be found in 
ref. [46]: Tounsi and Redlich show that the suppression factor governing the yields 
of strangeness 1, 2, 3 hadrons in proceeding from the canonical case (elementary 
hadron collisions and e + e~ annihilation to hadrons) to the GC limit (central AA 
collisions) is given by 

F = I S {V)/I Q {V) (5) 

where the / are modified Bessel functions of ascending order in strangeness S = 
1, 2, 3 and V is the quantum number coherent fireball volume that feeds downward 
to the frozen-out hadronic gas volume (which may initially be of equal size to the 
decaying fireball volume). The Bessel functions essentially depend on the coherent 
volume; the asymptotic limits for this ratio are 1 for V — > oo, and V/2 for V 
near zero. The ascend to the "no suppression limit" at large volume is steepest for 
singly strange hadrons, and slower for multiply strange species. At small volumes, 
the factor grows about linearly with the reaction volume. Fig. 7 refers to the top 
SPS energy, E/A=158 GeV corresponding to y/s=17.3 GeV per participant nucleon 
pair. From this calculation we infer an enhancement factor of two for singly strange 
hadrons in the large coherent volume limit, and a much larger factor for s=2 and 3 
hyperons. This study also shows that the canonical suppression - or grand canonical 
enhancement factor depends, in more detail, on y/s: it falls down with \fs increasing 
from SPS to RHIC energies. As we have seen this analysis captures the observations 
made at top SPS energy [30, 35]. However it also holds at lower SPS energy [41] and 
further downward to AGS [47] and even SIS [48] energy. Strangeness enhancement 
(GC behaviour) is, thus, not indicating deconfinement, as proposed by its pioneers 
[16]. There must be a more general property of the maximum density fireball stage 
causing GC behaviour: probably large size and coherence. On the other hand let us 
avoid a misunderstanding: a QCD quark-gluon plasma represents one such case of 
large size and quantum coherence, and if it decouples directly to a hadron gas one 
expects to observe GC order [12, 33]. But there ought to be other forms of fireball 
matter that decouple similarly. 

5.1 Fireball Size and Coherence 

In order to pin down the origin of GC freeze-out let us first look at the size depen- 
dence. At top SPS energy we have data for central C+C, Si+Si, S+S and Pb+Pb 
collisions from NA35 and NA49 [15,17,51]. Fig. 9 shows compilations of these data 
[51] concerning the average total ir + , K + , K~ and multiplicities of such central 
collisions (for Pb+Pb we include the top two centrality windows) . The ratio of kaon 
and <I> multiplicities to those observed for pions is plotted against the projectile nu- 
cleon participant number. We see that all these ratios shoot up steeply from the 
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(canonical) pp values that are included for reference [18]. The overall trend resem- 
bles the Tounsi-Redlich model [46] illustrated in Fig 7. We see that S+S already 
exhausts the GC limit to about 80%, with < N part >ss 55. In general the model 
predicts a transition from C to GC order that is yet steeper than the data. However 
we have to be precise here. Firstly, the model is run with constant temperature and 
baryochemical potential - perhaps not exactly true in the data. Second, < N par t > 
has a different meaning in the data and the model: the statistical model has in 
reality no such < N par t > parameter whatsoever, but only a reaction volume in 
which to enforce either local quantum number conservation or conservation on the 
average only (from the C to the GC situation). I.e. the authors of Fig. 7 [46] fit 
a coherence volume V, starting from pp analysis where they find V ~ 7/m 3 in a 
situation with, undoubtedly, two participants. With this normalization they merely 
relabel their volume by N par t = 2V/7 fm?. Comparison with AA data is thus not 
straight forward. We may read off Fig. 7 that for s=l the grand canonical limit is 
reached for "N par t" > 20 which thus really means V > 70/m 3 . Not unreasonable: 
the rms radius of 32 S is about 3/m, thus its area 19/m 2 . The rapidity distribution 
of K + has a FWHM of about 2.5 units [40]. Employing the Bjorken-equivalence 
[10] which equates the unit rapidity width with a spatial cylinder length of lfm we 
estimate the primordial source volume of a central S+S collision to amount to about 
50/m 3 . 
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Figure 9: The ratios of K + ,K~ and $ to positive pion production in full phase 
space, plotted vs. participant nucleon number, for central C+C, Si+Si, S+S and 
Pb+Pb collisions at 158 GeV per nucleon [15,17,54]. The values observed for p + p 
collisions [18] are given for reference. 



As a last, important step in our attempt to understand the onset and origin of 
grand canonical behaviour we shall demonstrate now that it is neither the partici- 
pant number nor the mere volume of the fireball by itself that matters (recall Fig.2). 
To this end we show in Fig. 10 (left column) the data of Fig. 9 confronted with the 
results of minimum bias PbPb collisions at similar energy. They exhibit a strikingly 
different, smooth ascend with N par t which thus can not be the right scaling with 
variable. 

We are now prepared to formulate the essential hypothesis that we have re- 
peatedly hinted at in sections 2-4. We propose that the crucial GC coherence ef- 
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feet resides in the space-time density of "successive" collisions at the microscopic 
level, during fireball formation time. This dynamical quantity is, of course, not an 
ingredient of the statistical model description of hadronization. But it might be 
the property that justifies a grand canonical description of the system right after 
hadronization. In order to obtain an estimate of the space-time density of "succes- 
sive" collisions in the early fireball, we have to consult a model of the microscopic 
dynamics. The calculation employs the Frankfurt UrQMD model [52]. 
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Figure 10: Minimum bias Pb+Pb data confronted with the data shown in Fig. 9 
(left column). Right: all data plotted versus collision density, ref. 51. Middle: 
the fraction of nucleons in each collision system undergoing more than one inelastic 
collision. 



These considerations are illustrated [51] in the right column of Fig. 10. The data 
of Fig. 10 (left column) are replotted here as a function of the average density of 
successive inelastic collisions that would characterize each of the collisions systems 
if a quasi classical cascade picture would be applicable. Important: we shall conclude 
shortly that this picture of sequentiality should not be applicable but we consider 
the implicit transport ansatz nevertheless to be reliable insofar as it captures the 
consequences of initial impact geometry and density distributions in a symbolic 
quantity "inelastic collisions per cubic fermi and unit time fm/c", in the emerging 
fireball. Fig. 10 shows a smooth increase of the strangeness enhancement quantities, 
with saturation occuring in central Pb+Pb. All colliding systems line up smoothly, 
and we see that saturation, i.e. the GC limit, is reached at < g >> 6. 

Now let us, finally, turn the argument around. A system with 6 separate inelastic 
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collisions concurrent, on average, in each of its space-time unit cells is a strange idea! 
The uncertainty relation tells us that the energy uncertainty at each "successive" 
step would exceed 1.2 GeV here. I.e. the average density of energy uncertainty 
would exceed the critical QCD energy density. Such a system can not possibly be 
described as a sequence of hadron collisions. Its decay must occur under global 
quantum mechanical coherence, from interfering local subprocesses (for which we 
have no microscopic model as of yet). This mechanism should be the origin of 
the apparent grand canonical coherence. The cascade-transport model itself thus 
indicates the point at which it becomes invalid (indeed it does not predict the GC 
multiplicity pattern). The new physics seems to set in with C+C already, and is 
well established from S+S onward (Fig. 9) where the hot core of the fireball may 
have a volume of about 50 fm 3 . Not terribly far from the prediction [46] of Fig. 7 for 
strangeness one. 

It would be extremely interesting to check also the "onset curves" for s=2,3 
with hyperon data for lighter collision systems. Meanwhile we conclude that the 
occurence of a global chemical potential in the grand canonical description reflects 
quantum mechanical coherence over a significant volume at the stage of hadroniza- 
tion decay of the fireball. It is important to repeat that in this view strangeness en- 
hancement is not, by itself, an unambiguous "smoking gun" signal of deconfinement 
and quark gluon plasma formation. High density hadronic matter (whatever that 
might be) should also exhibit grand canonical hadron yields. In both cases, however, 
strangeness enhancement signals the production of hitherto unknown QCD states 
of matter. This should even apply to the relatively modest effects observed in pA 
collisions [26, 27, 28] that we mentioned in sect. 2. The "tube" of target matter 
struck by the incident nucleon should also be a coherent (albeit small) system. 
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